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ABSTRACT 

An exploratory analysis has been made of the aeroelastic stability of the Space 
Shuttle Launch Configuration. The purpose of the analysis is to define critical flow 
phenomena with adverse aeroelastic effects and to develop simple analytic means of 
describing the time-dependent flow-interference effects so that they can be incoipo- 
rated into a computer program to predict the aeroelastic stability of all free-free 
modes of the shuttle launch configuration. 

Three critical flow phenomena have been identified: 

1, Discontinuous jump of orbiter wing shock 

2, Inlet flow between orbiter and booster 

3, H, O. tank base flow 

All involve hi^y nonlineai* and often discontinuous aerodynamics which cause limit 
cycle oscillations of certain critical modes. Given the appropriate static data, the 
dynamic effects of the wing shock jump and the HO tanlc bulbous base effect can be 
analyzed using the developed quasi-steady techniques. However, further analytic and 
experimental efforts are required before the dynamic effects of the inlet flow phenom- 
enon can be predicted for the shuttle launch configuration. 
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Section 1 
TNTRODUCnON 


The present analysis is the first step toward the development of a technicpie for 
the prediction of the aeroelastic stability of the Space Shuttle Launch Vehicle, The 
quasi-steady approach is similar to that used so successfully for the Apollo-Satum 
booster (Refs. 1-4). The technique was later mechanized (Ref. 5) and used routinely 
to predict the aerodynamic damping of the low frequency free-free bending modes of 
the various operational vehicles (Refs. 6-10), and was further extended to predict the 
gust penetration response (Ref. 11). The present analysis is speciSeally aimed at 
defining the critical, time-lag sensitive, flow phenomena and develop means of pre- 
dicting their dynamic effects. Areas in need of further analytic and/or experimental 
investigations are pointed out. 

Much valuable data are available in the literature (Refs , 12-15) on the rigid body 
dyna.mie stability of both orbiter and launch configurations . These data were invalu- 
able to the present analysis. The approach taken was to first examine the rigid body 
dynamic stability data for the characteristic signature of a flow field time lag; i.e, , 
opposite static and dynamic nonlinearities (statically stabilizing but dynamically un- 
damping or vice-versa). Next static force, pressure, and flow visualization data were 
analyzed to determine the causes of these dynamic anomalies. The phenomena were 
then modeled using quasi-steady techniques to predict the dynamic measurements 
from static wind tunnel data. 

Once a satistectory match with che rigid body dynamic data was obtained, the 
analysis was applied to the elastic vehicle. For each phenomenon the aerodynamic 
damping was computed for modes where it had an undamping effect. Tins gav(3 a 
representative indication of the seriousness of the various unsteady flow effects. 
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Because of the large number of modes that exist for the shuttle launch configuration 
only a few representative modes were analyzed. It cannot be known if these are worst 
case modes without analyzing all modes. However, a complete modal analysis will 
have to wait until the procedure is computerized for rapid calculation of a large num- 
ber of cases . 
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Section 2 

ORBITEE UNSTEADY AERODYNAMICS 


Hie orbiter wing is tbe major lift producing element of the launch configuration. 
Thus, it will dominate the elastic vehicle dynamics. It is natural, therefore, to first 
consider the imsteady aerodynamics of the orbiter alone without the complication of 
booster interference. Consequently, the analysis begins with the orbiter alone vehicle 
dynamics . 

2.1 Orbiter Aerodynamics 

Typical orbiter dynamic stability data (Refs . 12 and 13) are presented in Fig- 
ure 1. The opposition between static and dynamic derivatives is indicative of a Sow 
field time lag (Refs. 1, 2, and 16). The sharp peaks in the stability derivatives 
suggest the possibSity of hi^ly nonlinear or discontinuous static stability character- 
istics (Refs. 17 and 18). Oilflow photographs (Ref. 19) indicate three Sow types on 
the orbiter wing at transonic speeds (Fig., 2) , At low angler, of attack, the Sow is 
essentially attached with a very small Sow separation at the foot of the normal shock 
on the outboard wing (Fig, 2a). The forv/ard shock appears to emanate from the for- 
ward edge of a comer separation at the wing fuselage juncture whereas the QMS pods 
fix the location of the aft shock (Figures . 2a and 3a) . At some intermediate angle 
of attack the corner separation jumps to the strake apex, and the forward shock also 
jumps forward quite near the wing leading edge. The wing shock induced pressure 
rise couples with the leading edge suction to produce a large adverse pressure grad- 
ient which causes a large region of flow separation to occur over the outer wing panel 


^ 2 

Stable C + C • excursions correlate with unstable C - k C • excursions 

mq mcK m.Oi mq 

and vice versa. 

The QMS pods also effect the position of the foi*ward shock via the comer separation 
as will be discussed later. 
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(Figs. 2b and 3b). Finally, at a bigti angle of attack, the shock moves all the way to 
the leading edge and the entire outer panel Is stalled (Figs. 2c and 3c), 

It is evident that the flow at the wing-fuselage juncture somehow alters the con- 
ventional two-dimensional shock-boundary layer interaction. Lock and Rodgers 
(Ref. 20) have described the mechanism whereby conditions at the wing fuselage junc- 
tm’e determine the position of the forward wing shock. The flow component normal 
to the leading edge is accelerated by expansion over the upper ving and, as a result, 
the resultant flow is turned toward the fuselage. The flow must eventually be turned 
parallel to the fuselage. This is accomplished via a compression that coalesces into 
a single shock; the forward wing shock. The similarily to the orhiter wing flow is 
illustrated in Figure 4. It is the comer separation that furnishes the boundary condi- 
tion that fixes the forward wing shock on the orbiter. When the separation jumps to 
the stralce apex the boundary condition for the forward vdng shock changes suddenly 
causing the shock to jump too. 

The comer separation has a duality of character; (1) it acts as the boundary to 
which the leeside wing flow must adjust Itself, (2) it is vented through a vortex which 
generates Tift over the aft portion of the wing through vortex induced suction. Thus, 
it plays a fundamental role in producing all cf the orbiter nonlinear stability charac- 
teristics illustrated in Figures 5, 6, and T, 

The discontinuous shock movement caused by the comer separation jump gener- 
ates *the jumps in C and C at M = 1.2. Although the density of the data points 

nx 

is generally insufficient to aefine the nonlinearities, the normal force data strongly 
suggest a highly nonlinear or discontinuous behavior between OJ - 8“ and oj = 10°. 

A corresponding nonlinearity in (o!) is computed from the dynamic data (Refs. 12 
and 13), as will be discussed later. Following the jump remains stabiliaing 

and, as the vortex venting the comer separation grows with angle of attack, 
becomes more stabilizing while increases due to the vortex-induced suction 

over the aft wing (Refs, 21 and 22). Concurrently, plateaus as the vortex 

:K ■ 

Although the aft shocks have different origins . 
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suction opposes th© usual reduction of witli a. Finally, atMgh ct the wing 
stalls causing an unstable increnxent in C^, a reduction in and a decrease in 
C . , These trends are generally pervasive despite changes in OMS pod configuration 
(Figs, 5, 6, and 7 and Kefs, 23, 24, and 25, respectively). At M = 0,6, the data 
do not exhibit any discontinuity since no wing shock occurs at so low a Mach number. 
However, at a rj 8° an increase in and a more stable occurs dtje to 

the effects of the wing-fiiselage vortex. At high angles of attack, the stall nonlinearity 
is also evident. 


The center of pressure of each of these nonlinear incremental variations in the 
static stability data have been computed; i,e, , 



C — C 
mo* mao 

~ '^Ng'o 



where the subscript o denotes the slope at o' = 0, The results are shown in Fig- 
ure 8 . These data indicate that the nonlinear effects occur on the aft wing panel for 
both static and dynamic measurements. Both static and dynamic nonlinearities cor- 
relate with the change in fl'>w patterns from oilflow results (Fig, 9) , further verify- 
ing that the wing flow is responsible for both static and dynamic nonlinearities , 


The variation in the flow boundaries in Figure 9 are largely the result of varia- 
tions in OlVE pod configuration. Figure 10 shows that the wing flow patterns are 
sensitive to the OMS pod configuration (Kefs, 19 and C6). The OMS pods directly 
effect the position of the aft wing shock and indirectly they also effect the position of 
the forward shock. The OMS pods create an adverse pressure gradient at the wing- 
fuselage juncture due to the OMS ow dhock at high trsinsonic speeds (Fig, 11a) and 
due to stagnation at the OMS ii>.duced downwash at lower speeds (Fig, 11b), This 
feeds upstream through the thick viscous layer at the wing-fuselage juncture, affecting 
the location of the comer separation and thus the foarward wing shock. Unfortunately, 
the OMS pod coiifiguration is the most changed configuration feature of the orbiter. 

It seems to be different for nearly every wind tunnel test. Of the static tests, the 

The oilflow photographs were taken at a = 5° incrementsjthus the S'’ data bars in 
. Figure 9. 
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OMS pod configuration in Hef , 23 is the one most like the pods of the dynamic model 
(Refs. 12 and 13). 


2,2 Orbiter Rigid Body Bynamics 


All evidence indicates that the straite-fuselage vortex dominates the orbiter wing 
loads for a large angle of attack range , Whether it is the discontinuous shock jump or 
the suction effect of the vortex, the nonlinear wing load is a function of the corner 
separation. For a first approximation it was assumed that the crossflow at the strake 
apex determines the boundary layer condition at the wing-fuselage juncture and thus 
the corner separation extent. When the separation occurs at the strake apex, the 
situation is similar to that for free -body vortices on slender bodies of revolution 
(Ref. 27). For the shuttle this means that the crossflow at the strake apex determines 
the position of the comer vortex relative to the aft wing, thus determining the aft wing 
lift. At high angles of attack the strake and main wing vortices combine into one vor- 
tex (Fig. 3c). This rolled up vortex is sv/opt to the outboard wing panel where it bursts 
at stall. Thus, all the nonlinearities in the stability derivatives are determined by 
crossflow at the strake apex. Bi the dynamic case the induced load, the nonlinear 
increment, will lag the vehicle attitude due to the finite convection speed of the vor- 
tices (Ref, 28), That is, the induced pitching moment is a function of the angle of 
attack at the stralce apex at an earlier time (Refs, 6 and 22); 




m 




< 2 ) 


At 




U 


where x. is the apex station, x is the induced load station, and U = U is the 

A S 

convection speed of the vortex (Ref. 22). With the definitions in Fig, 12 


Cl = OtQ + 
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and 


(t-At) - e (t-At) - c ^ (t-At) 

0{t-At) = = 0 (t) cos(coAt)- i0 (t) sin(wAt) 

0 “ £j [i 0 (t)j 


therefore 

a 

0 (t-At) = 0 (t) cos(o)At)- Gm(wAt] 
and 

0 {t-At) “ 0 (t) cos (o>At) - sin (coAt) 
The aerodynamic inertia (0) terms are negligible. 


(4) 

(5) 

( 6 ) 


Combining equations (2), (3), (4) and (5) gives 


where 
and 

The induced damping derivative is added to the attached flow damping derivative 
for angles of attack where the strake-fuselage vortex occurs, L. e. , for a^CL 



a A^C 


m 




= - 2 AC 


mo; 


sin(coAt)^ / 

—L — cos(wAt) 


(7) 


w = 


CJC 

U 


A^C 


= C - C 
m m 
a OL ao 


( 8 ) 
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and for a < a 

V 




where 



2 C 


m 


ao 



( 10 ) 


( 11 ) 


is the trailing edge of the equivalent orbiter slender delta wing (Kefs, 21 and 22). 
Eqs. (O' and (10) can be used to estimate the orbiter damping as long as the induced 
pitching moment derivative is finite and constant; however, when the comer separation 
jumps, = w . An equivalent static moment derivative can be defined by inte- 

grating the discontinuous moment curve over the oscillation amplitude (Ref. 29). 



1 

A 9^ 



d 9 



+ 



A6 



A $ 


1 


( 12 ) 


provided that + A9 and - < gj^ - Ad 

where is the moment discontinuity and A6 the oscillation amplitude; likewise, 

the equivalent damping derivative is defined by integrating over the cycle (Ref, 29), 


ac- 

m 

^( 1 #) 


+ 27T 


7TA9 <j3 


COS Ipdip 




(13) 


where - wt and 2tt = coT; T is the period of oscillation (Fig. 13). 


This gives (following file method of Ref. 29); 



4AC , 


( 14 ) 
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where 


i/' = sin 


OJ-pv- OL 

D o 
A0 


cos j3 


and 


tan j9 = 2 


sin wAt , t 

— ?? — + cos coAt 


w 


CtJC 

2U 


The 3£nagnxtude of the moment discontinuity (AC ) was estimated from the spikes 
2 

in C - k C • from the dynamic test data (e,g. , at oj - 2 and 8 degrees in 

UGlO* Tyi(^ ^ 

Figures la and lb, respectively) using Eq. (12) (where A$ = 1®). These values 
were then substituted into Eq, (14) to obtain the damping spikes. The rest of the 
damping curve was computed using the static data of Eef , 23 since that OMS configu- 
ration best approjdraates the OMS configuration of the dynamic model. In Figure 13 
the results of these estimates are compared to the Langley damping measurements 
(Eef. 12). The agreement is gratifying since it verifies the flow model. Actually, in 
the stall region the time lag should probably be somewhat longer than that used in the 
present analysis as indicated by dynamic experiments (Eef. 30), This would give 
better agreement with tlie experimental damping values. However, stall occurs at too 
high an angle of attack to be of concern for the launch configuration. Thus, it is un- 
necessary to pursue the stall analysis further. 


*This is how the discontinuities in Figures 5-7 were obtained. 
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Sections 

lAUNCH GONPIGUEATTON DYNAMICS 


Armed with a satisfactory model of the unsteady aerodynamics of the orbiter wing 
one can now set out to analyze the dynamics of the launch configuration. 

3,1 Launch Configiiration Aerodynamics 

Rigid body pitch damping data for the launch confiLguration (Ref, 14) also exhibit 
the characteristic signature of a flow field time lag (Figs. 14 and IS). At M.- 1,2 
the large undamping peak at o; = 2. 2“ suggest a pitching moment discontinuity. 
Shadowgraph photographs indicate a region of flow separation on the HO tank due 
to the detached orbiter bow shock (Fig, 16), 

Pressure distribution data (Ref. 31) indicate how this separation responds to 
an^e of attack (Fig, 17). Crossflow over the HO tank nose forward of separation 
results in a thickening of the leeside boundary layer which in turn causes a forward 
movement of the separation point. The separation point location is a nonlinear func- 
tion of angle of attack (Fig. 18). The terminal normal shock standing between the 
orbiter and HO tank exhibits a nonlinear aft movement that correlates with the for- 
ward movement of the orbiter bow shock (Figs. 17 and 18), likewise, the static 
pressures aft of the shock on both the HO tanlc top and the orbiter bottom also 
exhibit a nonlinearity that correlates with the shock movement. Thus, the entire 
flow field is coupled. This is the first indication that the flow between the stages is 
analogous to inlet flow (Ref, 32) , 


In fairing the shock position curve the discontinuity has been positioned at a; - 2, 2 “ 
to correspond with the dynamic data. 
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The position of the normal shock between the stages is determined by the mass 
flow that can be accommodated by the more or less two-dimensional channel between 
the stages from the strake apex aft to the HO tank base. Any excess mass flow 
must be spilled laterally upstream of the stralce. The flow upstream of the orbiter 
bow shock is aligned axially and spilling begins aft of the shock. The spilling is a 
function of the pressure distribution over the top of the HO tanlc. When the pres- 
sures are high a relatively large amount of mass flow is forced out laterally. How- 
ever, when the pressure drops, as it does under the orbiter nose, the outflow 
decreases. When the orbiter bow shock is in its aft position the area of high pressure 
is reduced and less mass flow is spilled. In order to get rid of the excess mass flow 
the normal shock between stages moves forward, giving an additional high pressure 
region forward of the strakes, which spills the excess mass flow. Of course this does 
not happen instantaneously; a considerable time lag is involved. The effect of an 
angle of attack change must first be convected downstream to the orbiter bow shock 
before the separation point can move. The excess mass flow must then be convected 
into the "inlet'^ and the ”inlet" must fill before the normal shock is forced forward 
to the new equilibrium position. This filling time can be quite long, especially when 
one considers the considerable area for venting between the stages. 

3.2 Launch Configuration Rigid Body E^amics 

It is impossible to estimate the time lag theoretically. However, the combination 
of static pressure data, dynamic stability data, and quasi-steady theory supplies the 
tools necessary for estimating the time-lag. If both the static and dynamic effects of 
the discontinuity are known, Eq. (14) can be used to estimate the time lag. At. Tnis 
has been done for M = 1.2 using an iterative procedure. The jump increment was 
first talten in reference to the attached flow level to obtain a first estimate of the 
time lag. Then the plateau level for a > - 2 ° , on both sides of the undamping peak 
(Fig. 15), was computed. The jump was then incremented from the plateau and a sec- 
ond time lag determined. After just two iterations the results shown in Figures 19 


The attached flow damping was estimated using first order momentum theory (Ref, 33) 
for a slender body. Sfetic force measurements for the launch configuration were used 
to improve accuracy. 
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and 20 were obtained. The jump loads and 12ie effect of the continuous shock move- 
ment were estimated by integrating the pressure data of Ref. 31 (Figs, 21 and 22 
respectively). 

The very long time lag (the vehicle travels 26.2 or biter lengths during At) is 
not unreasonable for inlet flows (Ref. 32) , The good agreement with the Langley re- 
sults is strong evidence that the flow model is valid. 

The inlet flow effect is not restricted to M = 1. 2 but occurs over the Mach number 
range 0.9 < M ^ 1,4 (Figs. 23-25), However, at M = 0.9 the discontinuity in 
the shook position is not present since the shock stands so close to the shoulder. 

This explains why noundamping pealc appears in the M = 0.9 dynamic data (Fig. 14) 
although the plateau is present. The plateau is well predicted from the static induced 
derivatives (from the integrated pressure data of Ref. 31) with wAt = 60° (Fig. 26). 
Evidence of the strake-fuselage vortex effect is also present in the M = 0 . 9 data. 

The damping spdce is well predicted if one assumes that the wing shock jumps at 
Qt = 4.5®, rather than at oi = 2.2® as it did for the or biter alone. Evidently the HO 
tank and SRM's reduce the crossflow at the strake apex so the corner separation 
jump does not occur until a; =4.5®. Using the same ratio 

2,0 the shock jump will not occur untU. ci = 16° at M = 1,2 explaining why it was 
not observed in the booster damping data for oi :< 6° (Fig. 15 and Ref, 13). 

Another, simpler flow model was postulated earlier in the study. It was assumed 
that the induced loads did not involve an inlet-like flow. Rather, the orbiter bow shock 
was treated as a simple terminal normal shock with only boundary layer convection 
and accelerated flow time lag (Refe. 17 and 18), It was necessary to assume a large 
(.7 degree) amplitude hysteresis loop in order to predict the imdamping peak at 
M = 1,2 (Fig. 27). However, the plateau could not be accounted for, using the con- 
tinuous separation movement, because the time lag was too small. Although it is not 
conclusive, this certainly provides further evidence to support the inlet flow model. 
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Yaw damping data for the lamcli configuration with and without SKM's (Eef. 14) 
exhibit nonlinear damping characteristies (Fig, 28) that are indicative of bulbous base 
undamping with possible sting interference effecis (Refs, 34 and 35), A good approxi- 
mation to the yaw damping results is obtained for the configuration with SRM's off by 
simply carpet plotting the base increment fcom Ref. 36 to account for base radius and 
proximity to the sting flare‘'^,and adding the increment to the first order attached flow 
damping estnnate (Pig. 29). Of course the agreement for the configuration with SRM's 
on is not nearly as good since the SRM's certainly alter the base flow. 

Unfortunately, further experimental data are needed in order to separate the sting 
interference effect from the free wake bulbous base effect. Actually, three induced 
loads affect the dynamics of bulbous bases (Pig. 30), Porebody crossflow effects 
thicken the leeside boundary layer. This causes the lee side separation point on the 
base to move forward, creating a separation asymmetry that causes a statically de- 
stabilizing base load (Fig. 30a). This load is aerodynamically undamping due to the 
time lag associated with conveeting the boundary layer from nose to base. The free 
wake effect results from the pressure gradient across the inclined wake at 01 ^ 0 that 
eventually turns it in the streamwise direction. Ihe pressure gradient causes hi^er 
windward side wake neck pressures that are convected forward to the se'paration point, 
causing a separation asymmetry on the base that produces a statically stabilizing 
aerodynamically undamping base load (Pig, 30b), Finally, the asymmetric sting flare 
causes a wake flipping effect, that results in yet another statically stabilizing, undamp- 
ing, base load (Fig. 30c). 


The sting configuration used in r btaining Gie data in Refs. 14 and 36 is identical. 
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Section 4 

AEROELASTIC STABIIITY 


The equation of motion of the elastic vehicle describing single degree-of-freedom 
bending oscillations can be written (Refs ,, 1 and 2) . 


rsJ 

m 


q(t) + 2 CO -jSjj (D^ + D^) j q(t) 

+ co^ j^l . ^ (K^ + K^) 1 q(t)| = P(t) 


( 15 ) 


where P(t) is the buffeting force input. 


Ctoe requirement for stability is that the coefficient of q(t) cannot be negative. 


i.e. 


f <°S -0 <“> 

where D and D are the aerodynamic damping derivatives in separated and at- 

S 2 * 

tached flow respectively (a negative coefficient is damping). ^ is the structural 
damping as a fraction of critical and - puts the aerodynamic damping 

into the same form as the structural damping. For a discontinuous aerodynamic 
load the damping is obtained by integrating over one cycle similar to what was done 
for the rigid body damping (Refs, 17 and 18). 

27T+j3 

Cjj (l/>) cos (17) 


D - T =- 

s n’Aq Ct) 
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For the discoatinuous wing load illustrated in Figure 12 the instantaneous angle 
of attack is 

■ 

a^-ir e (18) 

The jump effect is the result of conditions at the stroke apex at the time t-At. 
Therefore 

(t-At) 

= ^o + 0^ (t-At) (19) 


with the elastic body coordinate system (Fig. 31). 

(fr-At) = <#)■ q (t-At) 

(t-At) = 4 (t-At) 


and 

q (t-At) - q cos (wAt) sin (ojAt) 

4 (t-At) = q cos (o?At) - ^ sin (wAt) 
where q = Aq sin n>t 


( 20 ) 


( 21 ) 


Therefore, 


or 


a. 

1 


= g;^ + (^^)Aq I sin (cot) cos (coAt) 

- cos(&jt)sin(coAt)j - ^ Aqco| cos (cot) 


Cl. - 


cos (co At) + CO sin (cot) sin (coAt)| 
+ sec 0 


( 22 ) 


(23) 


where 


^ = cot 
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and 




(3 =tan 


0' {^^) cos coAt- 0{|^) u 


(24) 


Thus , for the discontinuity that occurs at o:^ = ck^j 


= sin”^ 




Aq [0' (Ia^ cos &)At 0 (Ia^ U 


(25) 


The effect of the discontinuous orbiter wing load is found from applying Eq, (17) 
to the characteristics of Figure 12, i.e. 


(Tp) cos ipdip 



^s TTAq ct> 


+ I a (^) cos ^dip 


+ 


/ 


D 
27T+/3 


Cj^ {\p) cos 




where 


\ ° 


•^N = °N 

^2 ^^0!2 


% + ^ j 5 

% ■*■ ■*' “i > “d 


Integrating gives 

2 0(1 : c 

°s ^ ** "; ?Aqca ^ 


(25) 


(26) 


(27) 
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The aerodynamic damping is 



where 

°a =-% ‘ ^ (29) 

is the attached flow normal force derivative of the boost eoaSguratipn. 

Some discussion of the limits of applicability of the Q^uasi-steady theory is appro- 
priate before presenting results of the analysis. A simple description of the limits of 
tlie technique is given in Ref. 37 and Figure 32. It shows clearly that as long as the 
wavelength (A.) of the interference effect is large relative to the characteristic dimen- 
sion (e) of the effected body element the quasi-steady technique is valid (i.e. , when 
o/k .25). When the characteristic length approaches the i/2 wave length it is 
possible to get an upwash over the leading edge of the submerged body simultaneously 
with a downwash at the tail. Thus, the single lumped tail load approxima.tiou is no 
longer valid. Refining the approximation through the representation of the load dis- 
tribution fay several lumped loads is not the answer since the leading edge will alter 
the upwash (or wave shape) over downstream portions of the body. Thus, the appli- 
cability of the quasi-steady technique becomes questionable. Of cotirse the limit is 
hot rigid and the applicability of the results will gradually deteriorate as e/k 
exceeds 0,25. 


From Eq. (27) one can see that when the vehicle is describing infinitesimal am- 
plitude oscillations (Aq) at q;^ = with = Q the induced damping D^-+-oo. As 

the oscillation amplitude increases D_ becomes finite, limit cycle oscillations will 

s 

occur for modes such as the one shown in Figure 33 when the coefficient of q(t) hi 
Eq, (15) is zero. If the structural damping is zero, = 0 in Eqs. (16) and (28), limit 
cycle oscillations occur at D + D = 0 (Fig. 34). A more realistic value of the 
structural damping is ^ - .01; therefore, a limit cycle oscillation occurs where 


.01 


pus 

4com 


[°a + 


( 30 ) 
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Per tlie mode in Figure 33 limit cycle oscillations will occur for M = 0. 9 and 1. 2 
(Pig. 34). This limit cycle value is conservative in the sense that it may not have 
time to develop, but is unconservative in the sense that no coupling with other modes 
has been included. To determine exactly how large the maximum limit cycle will be 
requires a more complete modal analysis. 

Bather late in the study the work of Chipman and Bauch (Ref. 38) became available 
to the authors. Although their analysis revealed no aeroelastic instabilities at n = 0 
the location of the node line of the first torsional mode gives reason for concern in 
regard to the wing shock jump at a > 0. Figure 35 shows the nodal line superim- 
posed on M = 1.2 oil flow photographs. The shock moves forward of the node line 
with increasing angle of attack. The negative shock induced load will be aerodynamic- 
ally undamping when it is forward of the node (Fig. 26), The statically destabilizing 
jump load was damping for rigid body oscillations. However, the sign of the jump 
moment about the wing node is reversed and it will therefore cause undamping. It 
appears that the shock jump will be undampiag at M = 0.9 as well (Fig, 36), Time 
did not permit an analysis of the torsional mode; however it appears that at the very 
least limit cycle oscillations are possible in the range .9 < M < 1.2, when the wing 
shock jump occurs . 

The OIWK pod configuration has a strong influence on the wing shock position. 

Little or no flow separation occurs forward of the streamlined QMS pods; however, 
a very large separation occurs forward of the blunt QMS pods (Fig. 10) . That the 
QMS pods affect the position of the wing shock is shown in Figure 37. Evidently the 
vortices that vent the separation on the orbiter canopy energize the boundary layer 
on the top of the orbiter causing the separation to shrink and the orbiter wing shocks 
to move back with the separation. This is an effect of the blunt QMS pods that has 
not been analyzed as yet; however, the effect should be significant as it involves the 
orbiter wing load and the convection time of the vortices from canopy to OMS pod, 
which is nearly the same as the strake vortex lag. It appears that this effect can be 
analyzed with the same techniques used for the strake-fuselage effect on the wing 
shock. 
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The inlet lilce flow between orbiter and booster involves massive discontinuous 
loads which will certainly result in limit cycle oscillations for critical modes, Mg- 
ure 38 illnstrafes such a critical mode. Damping results are presented in 3?igL^re 39, 
Because of the very long time lag associated with the inlet flow even small errors in 
the time lag result in a large variatioa in the phase angle. Thus, for a 2.5% xmcer- 
tainty in the time lag (which is easily within the accuracy of the time lag estimates) 
the damping results cover the gambit from strongly damped to undamped with a limit 
cycle oscillation (Fig, 39), Unstable regions also repeat every 1/2 cycle (Fig, 40). 
The quasi-steady results are highly questionable in this case even though the technique 
is valid from the standpoint of interaction length to wave length ratios (c/\ < .25), 

The applicability is suspect because the ’’inlet walls” make 20 pulsations during the 
long time lag. Nevertheless, the possible implication of the inlet flow (possible struc- 
tural failure) is serious enough to warrant further investigation to assure structural 
integrity. 

The bulbous base effect has been analyzed for the mode shown in Figure 41. 
Unfortunately the free wake effect cannot be separated from the sting interference 
effect with the information currently available. Thus, the wake interference incre- 
ment was simply taken from the data of Figure 30 and applied directly to the booster 
mode (Figure 41) to obtain the results in Figure 42. Of course, the bulbous base ef- 
fect could be very sensitive to exhaust plumes, which have not been accounted for in 
these estimates . Judging by the amount of aerodynamic undamping predicted without 
plumes it seems prudent to investigate the bulbous base effect including the effect of 
exhaust plumes. This can be accomplished with currently available quasi-steady 
techniques given the static sting deflection derivative and plume increments. 


rf» 

The effects of exhaust plumes and control deflections must be included in such a 
study. 
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Section 5 
CONCLUSIONS 


An exploratory analysis of the space shuttle boost configuration has revealed three 
unsteady flow phenomena that cause aerodynamic undamping of certain critical free- 
free bending modes; 

1. TOng shock jump 

2. Met flow between stages 

3. HO tanic base flow 

Of the three effects items one and two involve discontinuous aerodynamics that 
result in limit cycle oscillations of certain critical free-free bending modes. The 
wing shock jump affects both booster bending modes and orbiter wing torsional modes. 

Limit-cycle oscillations threaten the structural integrity in two ways; either by 
outright over stressing the structure due to large modal deflections or via fatigue. 
Fatigue is a much more serious problem for the reusable shuttle orbiter than xt was 
for the single flight ApoUo-Satum vehicles. Since limit cycle oscillations will cer- 
tainly affect fatigue life a complete modal analysis is recommended to define the limit 
cycles for all critical modes. 

Given the proper static aerodynamic data the wing shock jump and the HO t an k 
base flow effects can be predicted using present quasi-steady techniques. However, 
the quasi-steady technique runs into accuracy problems for the inlet flow problem. 
Furthermore, the evidence supporting the inlet flow model is somewhat tenuous. It 
is, therefore, recemmended that the inlet flow be thoroughly investigated experiment- 
ally. Perhaps the most useful information could be obtained from a fluctuating pres- 
sure test of the flow between the stages . Coherence results could show conclusively 
the coupling of the inlet flow and would furnish accurate convection speed measurements 
that could be applied to the computation of the modal damping of the flight vehicle via 
quasi-steady techniques . 
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All three unsteady flow phenomena will be sensitve to Eeynolds number since they 
all involve flow separation. Since the quasi-steady theory relies on static wind tunnel 
data, which is usually obtained for Reynolds numbers two or three orders of magnitude 
below the flight value, the effect of Eeynolds number on the induced loads must be 
understood to allow extrapolations to Ml scale. It is recommended that the effects 
of Reynolds number on these critical unsteady flow mechanisms be investigated to 
determine proper scaling techniques . Additionally a means of simulating the fell 
scale Reynolds number in a subscale wind tunnel test by boundary layer trip design, 
location, etc, needs to be developed to properly simulate the fell scale unsteady flow 
effects. 
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App 0 ndix A 
NOMENCLATURE 


wing span 

reference length; mean aerodynamic chord 
damping derivative of elastic body, Eq, (15) 
spring constant, Eq. (15) 
free stream Mach number 

2 

pitching moment; coefficient, = m/(pU 2)/Se 

2 

normal force; coefficient, = N/(pU 2)/S 

2 

yawing moment; coefficient, C^=n/(pU 2)/Sb 
buffeting force 
pitch rate 

normalized coordinate 
yaw rate 

reference area; wing area 
time 

free stream velocity 
horizontal coordinate; (Fig. 12) 
vertical coordinate; (Fig. 12) 
angle of attack 
trim angle of attack 
yaw angle 
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A 

I 

K 

Q 

X 

P 

<P 

< 5 &' 

?/» 

CO, w 

A 

a 

D 

S 

T 

V 


increment 

dimensionless horiaontal coordinate, | = x/c 

structural damping, fraction of critical 

aerodynamic damping, fraction of critical 

pitch angle (Fig. 12) 

wavelength 

free stream density 

normalized modal deflection (Fig, 31) 

normalized modal slope (Fig. 31) 

cot 

circular frequency, co = coc/U 

Subscripts 
strake apex 
attached flow 
discontinuiiiy 
separated flow 

denotes trailing edge of equivalent delta wing 

vortex formation 

denotes zero angle of attack 

Superscripts 

induced e.g, = separation induced normal force 

d/dx 

d/dt 

d^/dt^ 
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Figure 1 Typical Orbiter Dynamic Stability Data 
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a) Q! = 5° 




c) O' = 25® 

Figure 2 Orbiter Leeside Flow Patterns 
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Figure 4 Comparison of Leeside Wing Flows 
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Figure 6 IVpical Static Data (Ref, 24) 
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Figure 8 induced Load Locations 
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Figure 9 Correlation of Data Nonlinearities with Oil Flow Results 
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Figure 10 Effect of OMS Pod CJonfiguration on Orbiter 
Flow Patterns at M = 1.2, a = 0 
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Figure 11 OMS Pod- Wing Flow fiiterference 





Figure 12 Defmition of Terms 
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Figure 14 Typical Launch Configuration Dynamic Stabiliiy Data at M = 0.9 
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Figure 16 Shadowgraph of Launch Configuration at M - 1.2 
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Figur© 17 Pressure Uistributions Between Sta.ges of Launch Configuration at M ~ 1,2 
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Figure 18 Correlation of Shock Locations and Pressure Jump at M = 1,2 
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Figure 20 Comparison of Qiiasl-steady Predictions with Laimch Configuration 
Dynamic Normal Force Derivative Measurements at M = 1,2 
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Figure 21 Jump Muced Load Distribution for inlet Flow, M = 1,2 
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Figure 22 Ihduced Norioal Force Derivative Distribution 
for Inlet Flow, M = 1 . 2 
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Figure 26 Effect of HO Tank Separation and Wing Shock Interaction on 
Lift-off Configuration ^gid Body Damping; M = 0 . 9 
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Figure 29 Effect of HO Tanic Bulbous Base on Launch Configuration D 3 mamic Stability 
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Figure 31 Coordinate System for Elastic Orbiter 
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Figure 32 limit of Quasi-Steady Theory 
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Figure 33 Sketch of Mode No. 14, f = 4,G5() Hz 
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Figure 34 Effect of Orbiter Wn.g Shock Jump on the Aerodynamic Damping 
of the Launch Configuration Mode No. 14 
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NODAL LINE 


Figure 35 Ck)r relation of Wing Torsional Mode Mode Line with 
Orbiter Wing Flow Patterns, M = 1.25 
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Figure 36 (Correlation of Wing Torsional Mode Node Line with Orbiter Mng 

Flow Patterns , M = 0 . 9, a = 0 


LMSC-D057194 






Figure 38 SIsetch of Mode No. 42, f = 9.120 Hz 
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Figure 41 SIcetch of Antisymmetric Mode No. 40, f = 8.43 Hz 
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^ = 8.^27 eps 



Figure 42 Effect of HO Tank Bulbous Base on Aerodynamic Damping 
of the Launch Configuration Mode No, 40 
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